Can one control systematic errors of QCD sum rule predictions for bound states? 
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We study the possibility to control systematic errors of the ground-state parameters obtained 
by Shifman-Vainshtein-Zakharov (SVZ) sum rules, making use of the harmonic-oscillator potential 
model as an example. In this case, one knows the exact solution for the polarization operator, which 
allows one to obtain both the OPE to any order and the parameters (masses and decay constants) 
of the bound states. We determine the parameters of the ground state making use of the standard 
procedures of the method of QCD sum rules, and compare the obtained results with the known 
£ — i exact values. We show that in the situation when the continuum contribution to the polarization 

, operator is not known and is modelled by an effective continuum, the method of sum rules does not 

f^*) ■ allow to control the systematic erros of the extracted ground-state parameters. 
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*0 ! 1. INTRODUCTION 

f4 ' A QCD sum-rule calculation of hadron parameters [l[ involves two steps: (i) one calculates the operator product 
Or expansion (OPE) series for a relevant correlator and obtains the sum rule which relates this OPE to the sum over 
O ■ hadronic states, and (ii) one extracts the parameters of the ground state by a numerical procedure. Each of these 
D ' steps leads to uncertainties in the final result. 

& • The first step lies fully within QCD and allows a rigorous treatment of the uncertainties: the correlator in QCD is 
not known precisely (because of uncertainties in quark masses, condensates, a s , radiative corrections, etc.), but the 
, corresponding errors in the correlator may be systematically controlled (at least in principle). 
^ ' The second step lies beyond QCD and is more cumbersome: even if several terms of the OPE for the correlator 
\ were known precisely, the hadronic parameters might be extracted by a sum rule only within some error, which may 
be treated as a systematic error of the method. It is useful to recall that a successful extraction of the hadronic 
parameters by a sum rule is not guaranteed: as noticed already in the classical papers [J , the method may work 
. . in some cases and fail in others; moreover, error estimates (in the mathematical sense) for the numbers obtained 
0^ ' by sum rules may not be easily provided — e.g., according to [1], any value obtained by varying the parameters in 
the sum-rule stability region has equal probability. However, for many applications of sum rules, especially in flavor 
physics, one needs rigorous error estimates of the theoretical results for comparing theoretical predictions with the 
experimental data. Systematic errors of the sum-rule results are usually estimated by varying the Borel parameter 
and the continuum threshold within some ranges and are believed to be under control. 

In this Letter we concentrate on the systematic errors of hadron parameters obtained by a typical SVZ sum rule. 
To this end, a quantum-mechanical harmonic-oscillator (HO) potential model is a perfect tool: in this model both 
the spectrum of bound states (masses and wave functions) and the exact correlator (and hence its OPE to any order) 
are known precisely. Therefore, one may apply the standard sum-rule machinery for extracting parameters of the 
ground state and test the accuracy of the extracted values by comparing with the known exact results. In this way the 
accuracy of the method may be probed. For a detailed discussion of many aspects of sum rules in quantum mechanics 
we refer to HUSH. 

2. THE MODEL 

To illustrate the essential features of the QCD calculation, let us consider a non-relativistic model with a potential 
containing both a Coulomb and a confining part, defined by the Hamiltonian 

H = H + V(r), H =?/2m, V(r) = a s /r + V con{ . (2.1) 

The polarization operator H(E) is defined by the full Green function G(E) = (H — E)^ 1 as follows: 

11(E) = (2ir/mf 2 (f f = 0\G(E)\n = 0). (2.2) 
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Fig. 1: The expansion of the polarization operator T1(E) in powers of the interaction: wavy lines (green) represent the Coulomb- 
part exchange, dashed lines (red) represent the confining-part exchange. 



The full Green function satisfies the Lippmann-Schwinger operator equation 

G- 1 (E)-G Q 1 (E)=V, 
which may be solved by constructing the expansion in powers of the interaction V: 

G(E) = G (E) - G (E)VG (E) + ■■■ . 



(2.3) 



(2.4) 



Respectively, for the polarization operator 11(E) one may construct the expansion shown in Fig. [T] If the confining 
potential is known, one can calculate each term of the expansion. If one does not know the confining potential precisely, 
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an explicit calculation is not possible. In this case one can explicitly calculate the contribution IIo + Ilg + LTg -\ 

of those diagrams which do not contain the confining potential (first line in Fig. [1]) and parametrize other diagrams, 
in which the confining potential appears, taking into account the symmetries of the theory. In QCD such diagrams 
are parametrized in terms of condensates and radiative corrections to them. 

In what follows we make a further simplification of the model and switch off the Coulomb potential, which we, in 
principle, know how to deal with. We shall retain only the confining potential, choose it as the HO potential 



V{r) = mLu 2 r 2 /2 



(2.5) 



and study H(fi) [the Borel transform of the polarization operator II(i?)], which gives the evolution operator in the 
imaginary time 



nO*) = (27r/m) 3/2 (f) = 0| exp(-H/fi)\fi = 0). 
For the HO potential (|2.5[) . the exact analytic expression for H(n) is known Q: 
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H(/i) 
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Expanding this expression in inverse powers of /z, we get the OPE series for n(/i): 

u? 19 ^ 
V + 480/7 



IIope(m) = Ho(m) + HiGu) + n 2 (/i) 
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(2.6) 



(2.7) 



(2.8) 



and higher coefficients may be obtained from (|2.7p . Each term of this expansion may be also calculated from 
and (12. 4|) . with n corresponding to Go- 

The "phenomenological" representation for II(^() is obtained by using the basis of hadron eigenstates of the model, 
namely, 



oo 

n(/i) = ^2 Rn exp(-£'„//i), 

n=0 



(2.9) 



3 



5(1 In) 



log noi) 




Rn 



exp(E ///) 




Fig. 2: Dashed (black) line - the calculation for the exact H(fi); dotted (red) line - the calculation making use of the OPE for 
II(/i) with 4 power corrections, solid (green) line - the same OPE with 100 power corrections. 



with E n the energy of the nth bound state and R n (the square of the leptonic decay constant of the n-th bound state) 
given by 

i?„ = (27r/m) 3 / 2 |*„(r = 0)| 2 . (2.10) 

For the lowest states, one finds from (|2.7|) 

E Q = ^cj, R = 2V2cj 3/2 , E x = 7 -lo, R x = 3\/2cj 3/2 , .... (2.11) 

Now, imagine we know n(/i) numerically and want to extract the parameters of the ground state E and Rq. At 
large values of the Euclidean time 1/p the polarization operator is dominated by the ground state. Therefore, one 
has a plateau in — gm^y logII(/i,) and in the product exp(I?o/ (Fig. [2]). In practice, working at sufficiently 

large values of 1/u, one constructs ~ -qhjja k>gll(//) and, after verifying that the plateau has already been reached, 
determines Eq. Then one studies II(u) exp(E /p) and reads off Ro. In principle, it is possible to reach the plateau 
if one retains sufficiently many terms of the truncated OPE: the radius of convergence of the truncated OPE to the 
exact function H(p) increases with the increase of the number of terms taken into account, and, e.g., for N > 100 the 
range of convergence already reaches the plateau (Fig. [5]) . 

In practice, however, one can calculate only the first few terms and, consequently, the radius of convergence is 
rather small. For instance, with 4 terms the convergence region extends only up to values of 1/ fi for which the ground 
state gives 85% of the correlator. 

The method of SVZ sum rules aims at extracting the parameters of the ground state from the first few terms of 
the OPE, e.g. in the region where the contribution of higher states is by far not negligible. 



3. SUM RULE 



The sum rule claims the equality of the correlator calculated in the "quark" basis and in the hadron basis 

uj 2 19 u 4 



Roexp(-E /n) + / rfzp phcn (z) exp(-z/u) = dzpo{z)exp(-z/(i) + u 3/2 
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4u 2 480 u 4 



, (3.12) 



with po(z) = 2y / z/n. Following 0], we use explicit expressions for the power corrections, but for the zero-order 
free-particle term we use its expression in terms of the spectral integral. The reason for this will become clear in a 
few lines. 

Let us introduce the effective continuum threshold z c g(fi), different from the physical /i-independent continuum 
threshold z cont , by the relation 



n coIrt (/-0 = / dzp phcn (z) exp(-z/u) 



dz po(z) exp(— z/fi). 



(3.13) 
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Generally speaking, the spectral densities p p hen{z) and po(z) are different functions, so the two sides of (]3 . 13[) may 
be equal to each other only if the effective continuum threshold depends on fi. In our model, we can calculate Il cont 
precisely, as the difference between the known exact correlator and the known ground-state contribution, and therefore 
we can obtain the function z e g(ii) by solving (|3.13[) . In the general case of a sum-rule analysis, the effective continuum 
threshold is not known and is one of the essential fitting parameters. 
Making use of (j3 . 13[> . we rewrite now the sum rule (|3. 12|) in the form 



R exp(-E /{i) = II(/i, Zoff (//)), 



where the cut correlator z e s(jj,)) reads 



(3.14) 



Zeff (P) 



H(n,Z e g(n)) = J dz^exp(~z/p) 



3/2 



19 oj 4 
480^4 



(3.15) 



As is obvious from (|3. 14[) , the cut correlator satisfies the equation 



d 



d(l/M) 



logn(/i,2 e ff(M)) = E . 



(3.16) 



The cut correlator is the actual quantity which governs the extraction of the ground-state parameters. 

The sum rule (|3.14[) allows us to restrict the structure of the effective continuum threshold 2 e ff (/•*)■ Let us expand 
both sides of Q3.14p near u>/fi = 0. The l.h.s. contains only even powers of y/ui/n; power corrections on the r.h.s. 
contain only odd powers of y/u/fi. In order that both sides match each other, the effective continuum threshold 
cannot be constant but should be a power series of the parameter y/ujjfj,: 



Zeff(p) 



UJ 

ZQ+ ZlJ — 



(3.17) 



Inserting this series in (|3.14[) and expanding the integral on the r.h.s., we obtain an infinite chain of equations 
emerging at different orders of y/uijfj,. These equations allow us to obtain for any Eq and Rq within a broad range of 
values a solution z e #(jj,, Rq, Eq) which exactly solves the sum rule (|3. 121) . Therefore, in a limited range of p, the OPE 
alone cannot say much about the ground-state parameters. What really matters is the continuum contribution, or, 
equivalently, z c s(p). Without constraints on the effective continuum threshold the results obtained from the OPE are 
not restrictive! 

The approximate extraction of Eq and Rq worked out in a limited range of values of p becomes possible only by 
constraining z c g(p). If these constraints are realistic and turn out to reproduce with a reasonable accuracy the exact 
Zcff(p), then the approximate procedure works well. If a good approximation is not found, the approximate procedure 
fails to reproduce the true value. Anyway, the accuracy of the extracted value is difficult to be kept under control. 
This conclusion is quite different from the results of QCD sum rules presented in the literature (see e.g. the review 
Q). We shall demonstrate that a typical sum-rule analysis contains additional explicit or implicit assumptions and 
criteria for extracting the parameters of the ground state. Whereas these assumptions may lead to good central values 
of hadron parameters, the accuracy of the extracted values cannot be controlled within the method of QCD sum rules. 



4. NUMERICAL ANALYSIS 



In practice, one knows only the first few terms of the OPE, so one must stay in a region of p bounded from below to 
guarantee that the truncated OPE series reproduces the exact correlator within a controlled accuracy. The "fiducial" 
[l[ range of p is the range where, on the one hand, the OPE reproduces the exact expression better than some given 
accuracy (e.g., within 0.5%) and, on the other hand, the ground state is expected to give a sizable contribution to the 
correlator. If we include the first three power corrections, III, II2, and II3, then the fiducial region lies at uj/ p < 1.2 
(see Fig. [3]). Since we know the ground-state parameters, we fix ui/p > 0.7, where the ground state gives more than 
60% of the full correlator. So our working range is 0.7 < u/p < 1.2. 

Obviously, if one knows the continuum contribution with a reasonable accuracy, one can extract the resonance 
parameters from the sum rule (|3. 1 2[) . We shall be interested, however, in the situation when the hadron continuum 
is not known, which is a typical situation in heavy- hadron physics and in studying properties of exotic hadrons. Can 
we still extract the ground-state parameters? 
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We shall seek the (approximate) solution to the equation 

oo 

Rexp(-E/n) + J dzp (z)exp(-z/fj,) = IIope(m) (4-18) 

in the range 0.7 < uj//j, < 1.2. Hereafter, we denote by E and R the values of the ground-state parameters as extracted 
from the sum rule (|4. 18|) . The notations E and Rq are reserved for the known exact values. 



A. ^-dependent effective continuum threshold. 



In many interesting cases the ground-state energy can be determined, for instance, from experiment. However, 
fixing the ground-state energy E equal to its known value Eg does not help: for any R within a broad range one 
can still find a solution z e ff(/x,i?) which solves the sum rule (|4. 18|) exactly (Fig. [3]). In order to extract R, one needs 
to constrain the effective continuum threshold; moreover, this constraint determines the actual value of R which one 
obtains from the sum rule. If one requires, e.g., z c ff (/i) > Eq for all then the sum rule (|4.18p may be solved for any 
R within the broad range 0.7 < R/Rq < 1.15 (Fig. [3]). Clearly, the sum rule alone, without knowing the continuum 
contribution, cannot determine R. 



B. Constant effective continuum threshold. 



Strictly speaking, a constant effective continuum threshold z c h (p) = z c — const is incompatible with the sum rule. 
Nevertheless, this Ansatz may work well, especially in our model: as can be seen from Fig. [3][a), the exact z e g(/x) is 
almost flat in the fiducial interval. Therefore, the HO model represents a very favorable situation for applying the 
QCD sum-rule machinery. 

Now, one needs to impose a criterion for fixing z c . A widely used procedure is the following Q: one calculates 



E(ji,z c ) = - 



d(l/fi) 



iog-n(/i, z c 



(4.19) 



which now depends on fi due to approximating z c g(p) by some constant. Then, one determines fio and z c as the 
solution to the system of equations 



d 

E(fio, z c ) = E , d^ E ^' z ^ 



0. 



(4.20) 



yielding z c = 2.454cj, fio/uj = 1 (Fig. Q|, and a good central-value estimate R/Ro = 0.96, with R being extremely 
stable in the fiducial range. 

Note, however, a dangerous point: (i) a perfect description of Tl(fJ.) with better than 1% accuracy, (ii) the deviation 
of E(n, z c ) from Eq at the level of only 1%, and (iii) extreme stability of i?(/x) in the full fiducial range leads to a 4% 
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Fig. 3: Different choices of the effective continuum threshold z e s(fj.) (a) and the corresponding R(fJ,) (b): 1 [solid (red) line] 
the exact effective continuum threshold as obtained by a numerical solution of (14.18|l for E = Eq and R = Ro. 2 [long-dashed 
(blue) line] the effective continuum threshold obtained by solving the sum rule (|4.18[) for E — Eq and R = 0.7Ro, 3 [dash-dotted 
(green) line] same as line 2, but for E = Eq and ii = 1.15-Ro- 
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E(ju)/Eo R(yu)/R 




Fig. 4: Constant effective continuum threshold z c : E(fj,) for three different values of z c (a) and the corresponding R(u) (b). 
Dotted (black) line - z c = 2.4 u, solid (red) line - z c = 2.454 u, dashed (blue) line - z a = 2.5 ui. 

error in the extracted value of R\ Clearly, this error could not be guessed on the basis of the other numbers obtained, 
and it would be wrong to estimate the error, e.g., from the range covered by R when varying the Borel parameter /j, 
within the fiducial interval. 



5. CONCLUSIONS 

Let us summarize the lessons we have learnt: 

1 . The knowledge of the correlator in a limited range of the Borel parameter /i is not sufficient for an extraction of 
the ground-state parameters with a controlled accuracy: Rather different models for the correlator in the form of a 
ground state plus an effective continuum lead to the same correlator. 

2. The procedure of fixing the effective continuum threshold by requiring the average mass calculated with the cut 
correlator to reproduce the known ground-state mass 0, Q is not restrictive: a /i-dependent effective continuum 
threshold z c s(fi), obtained as a solution of the sum rule (|4.L8|) . leads to the cut correlator (|3.15p which automatically 
(i) reproduces exactly E(fi) — Eq for all values of the Borel parameter /j,, and (ii) yields a /x-independent value of R 
which, however, may be rather far from the true value. 

3. Modeling the hadron continuum by a constant effective continuum threshold allows one to determine its value z c 
by, e.g., requiring the average energy E(n) to be close to Eq in the stability region. In the model under discussion 
(where the exact effective continuum threshold is, in fact, almost constant) one obtains in this way a good estimate 
R/ Rq = 0.96, with practically /i-independent R. The unpleasant feature is that the deviation of R from Rq turns out 
to be much larger than the variations of E(p) and R over the fiducial range. Therefore, the systematic errors cannot 
be estimated by looking at the variation of R within the fiducial domain of fi. 

Thus, we conclude that in those cases where the continuum contribution is not known (which is a typical situation 
encountered, e.g., in analyses of heavy-meson observables) the standard procedures adopted in a QCD sum-rule 
extraction of hadron parameters do not allow to control the systematic uncertainties. Consequently, no systematic 
errors for hadron parameters obtained with sum rules can be provided. Let us emphasize that the independence of 
the extracted values of hadron parameters from the Borel mass does not guarantee the extraction of their true values. 

Nevertheless, in the model under consideration the sum rules give good estimates for the parameter Rq. This may 
be due to the following specific features of the model: (i) a large gap between the ground state and the first excitation 
that contributes to the sum rule; (ii) an almost constant exact effective continuum threshold in a wide range of [i. 
Whether or not the same good accuracy may be achieved in QCD, where the features mentioned above are absent, is 
not obvious and requires more detailed investigations. One of the possible ways could be the application of different 
versions of sum rules to the same observable (see the discussion in [l(J [H| ) • 

Presently, this shortcoming — the impossibility to control the systematic errors — remains the weak feature of 
the method of sum rules and an obstacle for using the results from QCD sum rules for precision physics, such as 
electroweak physics. 1 



1 We would like to point out that with respect to the systematic errors of hadron parameters, the method of QCD sum rules faces similar 
problems as the application of approaches based on the constituent quark picture: for instance, the relativistic dispersion approach [T3l 
gave very successful predictions for form factors of exclusive D decays and provided many predictions for form factors of weak decays 
of B mesons |l3j . However, assigning rigorous errors to these predictions could not be done so far. 
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